Disorder and superconductivity : a new phase of bi-particle localized states 
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We study the two-dimensional, disordered, attractive Hubbard model by the projector quantum 
Monte Carlo method and Bogoliubov - de Gennes mean-field theory. Our results for the ground 
state show the appearance of a new phase with charge localization in the metallic regime of the 
non-interacting model. Contrary to the common lore, we demonstrate that mean-field theory fails 
to predict this phase and is unable to describe the correct physical picture in this regime. 

PACS numbers: 74.20.-z, 74.25-q, 74.40+k 
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The transition from superconductor to insulator (SIT) 
in the presence of disorder continues to actively attract 
the interest of experimentalists and theoreticians, alike. 
Experimentally, this transition has been observed in thin- 
films of various materials by several different groups, re- 
viewed for example in Theoretically, attempts were 
made to model the original fcrmionic problem by an effec- 
tive system of interacting bosons, investigated in depth 
both analytically and numerically While this ap- 

proach yields interesting results, it is important to study 
the fermionic models, which directly correspond to the 
experimental situation Q. Recently, fermion models 
have been studied in the context of the SIT by quantum 
Monte Carlo (QMC) methods |j| which reproduce the 
physics accurately. This method yielded the transition 
from superconducting to insulating behavior as a func- 
tion of disorder strength. However, several unresolved 
questions pertaining to the physical origin of this tran- 
sition remain open. Further, a more microscopic, qual- 
itative understanding of the effects of disorder and the 
properties of the localized phase would be desirable. 

For weak disorder strengths, Anderson's theorem || 
guarantees that the superconducting phase persists de- 
spite disorder and the superconducting gap and other 
thermodynamic properties remain unchanged. The case 
of stronger disorder and non-uniform order parameters 
can be treated within the Bogoliubov- de Gennes ap- 
proach Such studies have been carried out by sev- 
eral groups ||]§] an d the study for s-wave superconduc- 
tors shows the persistence of a spectral gap with rela- 
tively strong disorder Q. In the spirit of Anderson's 
theorem, it is expected that the superconducting phase 
penetrates the localized non-interacting phase until the 
disorder strength is sufficiently large to overcome the su- 
perconducting gap. In this framework, the boundary be- 
tween the two phases can be estimated from the relation 
A ~ Ai ~ (y F l d )~ x } where A is the BCS superconduct- 
ing gap, Ai is the level spacing inside a grain of one par- 
ticle localization length I in the localized non-interacting 
phase and vp is the density of states at the Fermi energy. 
However, it is not clear whether the mean- field approxi- 
mation (MFA) remains valid in the presence of medium 



to strong disorder and its validity is worth testing, by 
comparisons with QMC calculations. Indeed, an indica- 
tion of the failure of the MFA comes from recent studies 
of the Cooper problem for two interacting particles above 
a frozen Fermi sea, in the three dimensional (3D) An- 
derson model (incorporating disorder) Jl(|. This study 
showed that the attractive Hubbard interaction creates 
localized pairs in the non-interacting metallic phase, in 
which the one-particle eigenstates are delocalized. This is 
in contradiction to the predictions of the MFA according 
to which Cooper pairs are delocalized in this regime. This 
disagreement is expected to be even stronger in lower di- 
mensions (namely, the 2D case). Further, the Cooper 
approximation cannot capture all the many-body effects 
and this renders interesting the study of the full, inter- 
acting, disordered fermion model with finite density of 
particles. 

Contrary to the expectations arising from Anderson's 
theorem, we find that the attractive interaction that leads 
to superconductivity gives rise to a phase of bi-particle 
localized states (BLS) in the metallic regime of the non- 
interacting disordered model. This confirms the indica- 
tions obtained from the generalized Cooper problem |lC| ] 
discussed above. Thus, our results, obtained by exact 
treatment of all many-body quantum effects in a realistic 
model at finite particle density, demonstrate convincingly 
the existence of the BLS phase. Furthermore, we show 
that the Bogoliubov - de Gennes MFA results are qualita- 
tively incorrect in this phase. The BLS phase is observed 
at medium disorder strength, while in the limit of weak 
disorder, the system is superconducting, as expected. 

To investigate the interplay of interactions and disor- 
der in the SIT regime, we have studied the attractive 
Hubbard model with disorder, on a square lattice. The 
model Hamiltonian reads, 



H = H A + Hj 

= \- t E & i,aH<r + E ei4<A<r 

(ij),a i,a 



where the a\ a (&i t(T ) are the creation (annihilation) oper- 
ators for a fermion of spin a at site i with periodic bound- 
ary conditions, hi a is the number operator for spin a at 
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site i, t is the hopping parameter, the Hubbard parame- 
ter, U < 0, gives the strength of the screened attractive 
interaction and e,, the energy of site i is a random number 
drawn from a uniform distribution [— W/2, W/2], which 
parameterizes the disorder. The first two terms represent 
the Anderson Hamiltonian and the last term represents 
the interaction Hi. The filling factor v = N p /(2 x L 2 ), 
where N p is the number of fermions (particles) and L the 
linear dimension of the system and the total number of 
sites is L 2 . It is known that the attractive Hubbard inter- 
action induces superconductivity; thus, we have a useful 
model without inquiring into the physical origin of the 
pairing interaction. 




FIG. 1. Charge density difference, 5 pi (see text for def- 
inition), at sites i of Hamiltonian (1), for a 12 x 12 lattice 
with W/t — 5 and N p = 72 fermions for a single disorder 
realization. The upper part (a) is with U/t = and the lower 
part (b) is the result from a Bogoliubov-de Gennes calculation 
with U/t = -6.0. 

We have studied this model by the projector quantum 
Monte Carlo (PQMC) method, which has no fcrmion sign 
problem for U < 0. As is well known this method 



treats the many-body problem exactly up to statistical 
errors and gives direct access to the ground state proper- 
ties of the system. The systematic error arising from the 
discrete symmetric Trotter decomposition, of step At, is 
of order (At) 3 . The simulations were carried out with 
At = 0.1, with up to 60 time steps, in the S z — sector. 
The filling factor v was tuned around quarter filling, for 
different disorder and interaction strengths, varying the 
system size L. The largest system studied was 144 sites 
with 74 particles. We note that this model has been stud- 
ied in the grand canonical ensemble by a finite temper- 
ature QMC method on an 8 x 8 lattice with v = 0.4375 
[f|. The SIT was estimated in this work to occur at a 
critical value of disorder strength W c « 3.5t. However, 
the system size dependence of the observables was not 
investigated in their work. 




FIG. 2. Charge density difference, Spi at sites i of Hamil- 
tonian (1), obtained from QMC calculations, shown here for 
a lattice of 12 x 12 sites, with N p = 72 fermions, for the 
same disorder realization as in Fig. 1. The upper part (a) is 
with W/t — 2.0, U/t = —6.0 and the lower part (b) is with 
W/t = 5.0 and U/t = -6.0. 
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To understand the physics of the model, we consider 
the difference in charge density in the ground state, 5pi — 
Pi{N p + 2) — pi(N p ), where pi is the charge density at site 
i and N p — L 2 /2 corresponds to quarter filling 8pi = 
2). For U = 0, the charge density difference Spi is given 
by the one-particle probability fi — \ifji\ 2 of the eigenstate 
at the Fermi level (fi = Spi/2). Due to this relation, the 
analysis of this characteristic, even in the presence of 
interactions, allows us to determine whether the charge 
is localized or delocalized at the Fermi level. 

An example of this characteristic for a typical disor- 
der realization is presented in Fig. 1. The results clearly 
show that the eigenstate at the Fermi level is delocalized 
in the non-interacting case at the given disorder strength 
W/t = 5. Of course, in the limit of L — > oo, eigenstates 
are localized for U = Nevertheless, for the finite 
system sizes in Fig. 1, the localization length is larger 
than the system size and the eigenstates correspond to a 
metallic regime. Certainly, it would be desirable to study 
this problem in the 3D case, where the Anderson transi- 
tion clearly separates the non-interacting delocalized and 
localized phases. However, it is presently too expensive, 
numerically, to study comparable system sizes in 3D. In 
spite of this, our data show that important physical in- 
formation can be obtained from the 2D case. 

The results obtained from the Bogoliubov-de Gennes 
MFA for the same disorder realization, with interaction 
strength U/t = —6 are shown in Fig. lb. These calcu- 
lations were carried out as described in [|]||. The eigen- 
vectors and the quasiparticle excitation energies were ob- 
tained self-consistently, with good convergence. In addi- 
tion, we obtained the BCS value for the order-parameter 
A w 1.36* in the limit of W/t = 0, at L = 12 and 
U/t = —4, as in ||. In the limit U = 0, this method 
reproduces Spi for finite disorder strengths W/t. A com- 
parison of Figs, la and lb indicates strongly that within 
the MFA, the charge density at the Fermi level is highly 
delocalized by the introduction of interactions. This 
result indicates that interactions smoothen out charge 
fluctuations within the MFA. Thus, at W/t = 5 and 
U/t = —6, the MFA gives completely delocalized, metal- 
lic behaviour. 

The PQMC results for Spi for the same disorder re- 
alization as in Fig. 1 are displayed in Figs. 2a and 
2b, corresponding to two values of the disorder strength 
(W/t = 2 and 5). The model parameters used are identi- 
cal for Figs, lb and 2b. It is apparent from Fig. 2b that 
the charge at the Fermi energy is localized. Thus, the 
PQMC result, with a proper treatment of interactions 
and disorder, differs qualitatively from the Bogoliubov - 
de Gennes mean-field calculation, which gives the wrong 
physical picture. We also emphasise that for U/t = 0, 
the charge variation Spi is delocalized (Fig. la) and the 
pronounced, localized peak of charge in Fig. 2b is not 
correlated to the charge distribution of Fig. la. With 
decrease of disorder strength W, this localized charge 



peak seen in Fig. 2b disappears as shown in Fig. 2a. In 
fact, the charge distribution in Fig. 2a becomes closer 
to that of the corresponding delocalized, non-interacting 
case (not shown). These data show that the SIT takes 
place in the metallic non-interacting phase. 

For a more quantitative description of this unusual 
SIT, we have studied an effective inverse participation 
ratio (IPR), £ for an added pair of particles. The IPR 
is defined as £ = (X^(^/°i/2) 2 } -1 : where the brackets de- 
note the average over 16 disorder realizations. For U — 0, 
this definition exactly reproduces the one-particle IPR at 
the Fermi level. In the presence of interactions, this quan- 
tity allows us to study the charge distribution of added 
pairs at the Fermi level, beyond the Cooper approxima- 
tion used in jl(J. Physically, the IPR gives the number 
of sites visited by a pair. Thus, the evolution of the IPR 
with system size L, can be used to determine the transi- 
tion from delocalized to localized behavior. 




FIG. 3. PQMC result for the IPR £ obtained from the 
charge density distribution of the added pair vs. linear di- 
mension of system L. Dashed lines are for W/t = 2 and full 
lines for W/t — 5, with U/t = (circles), —4 (squares) and 
—6 (diamonds). The average is carried out over 16 disorder 
realizations. 

The PQMC results for the IPR as a function of sys- 
tem size L, disorder and interaction strengths are shown 
in Fig. 3. The data clearly show two main phenomena. 
Firstly, it is clear that interactions always diminish the 
IPR, compared to the non-interacting case. Secondly, 
for moderate disorder strengths the pair becomes local- 
ized. Indeed, for W/t = 5, £ continues to grow with L 
for U/t = (metallic non-interacting regime), while it 
remains constant for U/t = —6. For example, in the 
latter case, £ = 6.7 ± 0.8 is much smaller than the to- 
tal number of sites (L 2 = 144), clearly indicating pair 
localization induced by the attractive interaction. This 
is in great contrast to the MFA result obtained as de- 
scribed for Fig. lb, which gives complete derealization 
with £ w 124. In fact, the £ from the MFA is signifi- 
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cantly larger than the non-interacting value of £ = 33.0, 
that once more shows that the MFA gives a physically 
wrong answer. With the finite system sizes of our sim- 
ulations, it is difficult to determine the position of the 
SIT precisely. However, we note that for W/t = 2, the 
system remains delocalized even for U/t = —6. Thus, we 
estimate that for U/t ~ (—4,-6), the critical disorder 
strength is W/t ~ 4. This is in satisfactory agreement 
with the results obtained for fixed system size L — 8 in 
H . We note that, clearly, well below the SIT, at weak dis- 
order, the Bogoliubov - de Gennes theory becomes valid. 

In fact, the appearance of this BLS phase can be un- 
derstood from the following heuristic physical argument 
fiof . The attractive Hubbard interaction strongly favors 
pair formation. This effectively doubles the mass of the 
charge carriers m* oc l/t and thus enhances the effects 
of disorder induced localization, which are proportion to 
W/t oc m* . We believe this to be the correct physical 
reason for the appearance of this BLS phase. 

In conclusion, our studies based on the PQMC method 
show the appearance of a localized phase of pairs which 
appears in the metallic regime of the non-interacting sys- 
tem. The comparison of these results with Bogoliubov - 
de Gennes mean-field calculations show that the MFA 
does not predict the existence of this phase. As the 
non-interacting states are delocalized for the moderate 
disorder strengths corresponding to the BLS phase, ar- 
guments in the spirit of Anderson's theorem should def- 
initely drive the system superconducting at the expense 
of the non-interacting localized insulating phase. On the 
contrary, our results provide evidence for the appearance 
of a new phase which is unexpected from the accepted 
viewpoint. We expect that this effect exists as well in 
three dimensions where it should be much more pro- 
nounced due to the sharp Anderson transition between 
the non-interacting metallic and insulating phases. 
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